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Abstract 

A weak scattering potential imposed on a CuOi layer of a cuprate superconductor modulates 
the local density of states N(x,u>). In recently reported experimental studies Q, Q,Q], scanning- 
tunneling maps of N(x,u>) have been Fourier transformed to obtain a wave- vector power spectrum. 
Here, for the case of a weak scattering potential, we discuss the structure of this power spectrum and 
its relationship to the quasi-particle spectrum and the structure factor of the scattering potential. 
Examples of quasi-particle interferences in normal metals and s- and d-wave superconductors are 
discussed. 

PACS numbers: 74.25.Jb, 74.25.-q, 74.50. +r 
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I. INTRODUCTION 



A weak scattering potential imposed on the C11O2 layer of a cuprate superconductor cre- 
ates ripples in the local tunneling density of states N(x, u) due to quasi-particle interference 
scattering. It was suggested that scanning tunneling measurements of the spatial and fre- 

[uency structure of N(x, u) could provide information on the k and oj dependence of the gap 
|. Recently, the introduction of high-resolution Fourier-transform scanning tunneling mi- 
croscopy (FT-STM) has provided a powerful new technique for studying this. In this 
approach, a BISCO crystal is cleaved exposing a B1O2 layer. Then an STM measurement 
of the local tunneling conductance dI(V, x)/dV is taken over a predetermined L x L grid of 
points that cover a region of order 600Ax600A. Assuming that the tunneling conductance is 
proportional to the underlying density of states of the Cu02 layer these measurements 
give an STM map of the local tunneling density of states N(x, to) with uj = eV. This map 
is then Fourier transformed, 

N(q,u)= e- i ^N(x l ,u), (1) 

xie(LxL) 

and the wave-vector power spectrum, 

|A(g,^)| 2 

P(q,v) = — — 2 — , (2) 

determined. Typically, the square root of the power spectrum, which is proportional to the 
magnitude of N(q,u), is plotted and we will follow that practice as well. Here, we will 
discuss the structure of P(q,u) and its relationship to the quasi-particle spectrum and the 
structure factor of the scattering potential. 

For an isotropic s-wave superconductor with a circular normal state Fermi surface, the 
ripples in N(x,u) produced by a weak scattering center form a circular pattern whose 
amplitude and wave length depend upon the bias voltage u> = eV. However, if the gap has 
d x 2_ y 2 symmetry, the ripples emanating from a scattering center appear as a characteristic 
set of rays whose wave length and amplitude vary with their angular direction and the size 
of the bias voltage u = eV)^, |6(. Two examples of this are shown in Fig. 1. It is the STM 
measurements of these modulations in N(x, to) that provide information on the wave vector 
and frequency dependence of the gap. Indeed, in their FT-STM power spectrum study of 

n n 

BISCO, Hoffman et. al [1( and McElroy et. al |3[ found frequency-dependent structure in 



P(q,u) which they argued were consistent with the Fermi surface dependence of A(k) as 

q n n 

seen from ARPES measurements |7(. In this work |1|,|3|, these authors suggested that the 
FT-STM data could be analyzed in terms of a set of frequency-dependent wave vectors 
q a (u;) which connect the tips of the constant energy contours specified by the quasi-particle 
dispersion relation 



The contours of the solid regions shown in Figs. 1(b) and (d) are the constant quasi-particle 
energy contours for a cylindrical Fermi surface e& = k 2 /2m — fi with a gap A& = AOCOS26 1 
for u) = 0.5A and 1.1A , respectively. Also shown in Fig. 1(b) are several of the q Q (co>) wave 
vectors introduced in Ref. |3j. The wave vector q'(u;) is an additional wave vector that we 
will discuss. The wave vectors qi(o>) and qs(w), along with their symmetry- related counter 
parts (not shown), are associated with the structure of the ripples in N(x, u) seen along the 
x and y axes of Fig. 1(a). Likewise, the wave vectors Cfa(u)), q'(^), and cfr(uj) determine the 
structure of the ripples along the 45° axes of Fig. 1(a). In conjunction with this experimental 
work, Wang and Lee ^] recently reported numerical calculations for the case of a single 
impurity which clearly showed the quasi-particle interference arising from qi(o>) and q7(o>). 
In addition, as u is varied these calculations showed that a rich, kaleidescope-like, structure 
appears in the wave vector power spectrum when it is folded back into the first Brillouin 
zone. 

In a similar way, for u = l.lAo, the wave vectors Q a shown in Fig. 1(d), determine the 
structure of N(x, u) seen in Fig.l (c). Here the ripples along the x and y-axes are associated 
with Q2 and Q3 (and their symmetry related counterparts in the y direction), while those 
along the diagonal are associated with Q\ (and its counterparts). In practice, when uj > Aq, 
inelastic scatering leads to a damping of the ripples in N(x,u), making the structure in 
\N(x, ui)\ associated with the Q a wave vectors difficult to detect. 

Now, the quasi-particle interference pattern shown in Figs. 1(a) and (c) are for a single 
impurity. For a particular surface region over which the STM measurements are made, there 
will be an array of scatterers leading to a complex overlap of ripples. Here we will discuss 
how the Fourier transformed wave vector power spectrum P(q,u) allows one to disentangle 
the quasi-particle interference effects from the static structure factors of the scatterers. In 
Section II we show that for the case of a weak scattering potential, P(q,u) factors into 
one piece which contains information on the nesting properties of the Fermi surface times a 



piece which is proportional to the static structure factor of the scatterers. We also note that 
these measurements contain information on the one electron self-energy. Various examples 
are analyzed to show the type of information that is in principle contained in the FT-STM 
data. In Section III, the case of a layered 2D superconductor is studied with results for both 
s-wave and d-wawe gaps discussed. Section IV contains our conclusions and Appendices A 
and B contain more details of the calculations. 



II. FT-STM POWER SPECTRUM OF A 2D NORMAL METALLIC LAYER 

To begin, we first consider the case of a normal metallic 2D layer. Suppose it is exposed 
to a weak potential 

V = J2 [ d 2 x5e{x)i)\(x)i) s {x). (4) 

s J 

with 5e(x) an energy shift at x. For a BISCO-like system this local energy change 5e(x) 
could arise from secondary effects associated with disorder away from the C11O2 plane, the 
regular potential of the B1O2 layer or possibly a weak static stripe potential. This interaction 
creates a ripple in the one-electron Green's function leading to a modulation in the local 
tunneling density of states. For the case of a weak potential which we will focus on, a Born 
approximation is appropriate so that the single particle Green's function is given by 

G{x, x', to) = G (x- x', u) + J d 2 x" G {x - x", u) 5e(x") G Q (x" - x', u). (5) 

Here, G Q (x, u) is the Green's function of the unperturbed system. Then the change in the 
single spin tunneling density of states at position x is given by 

5N(x, u) = ~Im J d 2 x" G Q {x - x" , u) 5e (x") G {x" - x, u) (6) 

Taking the spatial Fourier transform of 5N(x,u), on the L x L grid of points {x^} specified 
by the STM measurements, one finds for q ^ 0, that 

N{q,u)= V e- iq - Xl 8N(x h u) = J m A(g,w) (7) 

xi^(LxL) 



with 



A ( g? u) = J £x e iq - x G (x, u) G (-x, u) (8) 
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and 

5e(q)= [ d L^6e(x)e- i ^. (9) 



a 2 



Here a is the lattice spacing of the L x L STM grid and q = 2n(n x , n y )/La with n x and 
n y integers running from —L/2 + 1 to L/2. The wave- vector power spectrum of the local 
tunneling density of states, Eq. (J2J), is therefore given by 



—Ira A (q, u) 

7T 



2 

2 



\5e(q)\\ (10) 



Here N = L x L is the number of sites in the sampled region. Thus, in the weak scattering 
Born approximation, P(q,u) separates into a piece \ImA (q,uj)/7r\ 2 which describes the 
quasi-particle interference and a piece |5e(g)| 2 which is the static structure factor, S(q), of 
the scattering potential. 

For the random impurity case, one could imagine making STM maps over a large number 
of different Lx L regions. Then by averaging the structure factor over these maps one would 
obtain 

«^e* (,*0) (11) 

where n, is the area impurity concentration and Se is a site energy shift. In this case, P(q, u) 
would simply be proportional to n^e 2 times the quasi-particle interference factor from a 
single impurity. However, this is not the way the experiments are done. Rather, a single 
STM map on a finite Lx L grid of points covering a specific region is measured. In this case, 
|5e(g)| 2 versus q exhibits fluctuations bluring the image of \ImA(q, u)/tt\, although one can 
still resolve structure in false color 2D (q x ,q y ) maps of \N(q,u)\. However, as discussed in 
Appendix A, by averaging the power spectrum over blocks of width (Aq x ,Aq y ) about each 
q, the fluctuations can be reduced if the impurities are randomly distributed. Naturally, 
this reduces the momentum resolution. However, if the change in q of the quasi-particle 
interference response is predominately along a given momentum direction, one can average 
over a region of q values perpendicular to the direction of interest, reducing the fluctuations 
but maintaining the q resolution in the direction of interest. Here we will assume that a 
suitable average has been done and use the impurity structure factor given by Eq. (fTTj). 
Appendix A contains a further discussion of the effect of impurity induced fluctuations. 

For the case in which the scattering occurs from a regular lattice such as the Bi lattice, 
one has 

Mq)\ 2 = 5e 2 N% tGn (12) 



Here G n is a reciprocal lattice vector of the Bi lattice along with the satellite wave vectors 
associated with the supermodulation of the BiOi layer. One could also have a "random" 
array of stripe domains with 

(|fcW|2> = N - L-gly + p + fa-g^ + p + ( ' h ^ - q " q « ^ A s ' (13) 

Here Q x = 2n/l x with l x the stripe spacing, 2n /V is the characteristic size of a domain and 
Ni is the average number of domains in an L x L region. Here we have taken only the 
first Q x harmonic. The form factor associated with the charge distribution of the stripes 
suppresses the response at higher multiples of Q x . We will examine the effect of an array of 
scattering centers in Section III. 

Turning next to the quasi-particle interference response, we begin by looking at G (x,u>) 
for a free 2D electron gas. In this case, for u > 0, 

7^ 6 T • r = -i vriV(O) (k(u)r) , (14) 
(2-Ky u - e k + id 

with 6^ = k 2 /2m — fi and mu = k 2 F /2m. Here N(0) = m/27T is the single spin electron 
density of states for the 2D free electron gas, Hq is the zeroth order Hankel function of the 
first kind, r = Iscl and 



k(u) = k F ^l + ^-. (15) 

When k(uj)r is large 

i 

G o (x,u;)~-iN(0) (^) 2 e<^>-%). (16) 

and the spatial modulation of 5N(x, u) which varies as the square of G D (x, u) is characterized 
by a wave vector 

q(u)=2k F ^l + ^. (17) 

The quasi-particle interference response function A(q, uj) for the 2D electron gas is calcu- 
lated in Appendix B. The result of this calculation shows |2j that 

, 8^(0) for > 2k{u) 

-ImA(q,u) = { V^-^H ( 18 ) 



7T 



for q < 2k(u) 



Thus, the wave-vector power spectrum for the 2D free electron gas has a cusp at q equals 
2k{u). ImA(q,u) vanishes for q < 2k(u) and diverges as (q — 2k(u))~ 1 / 2 as q approaches 



2k(u) from larger values. As noted in the appendix, ReA(q,u) has a similar cusp as q 
approaches 2k(u) from below. Basically, there is just a shift of phase of tt/2 in A(q,u) when 
q passes trough 2k(u). If the impurities are dilute, but the scattering from a given impurity 
is strong, one still has P(q,cj) proportional to the impurity concentration. However, in this 
case jlOj], because of the phase of the t-matrix one will have singularities on both sides of 
2k{u). 

The experimental FT-STM data has been reported as the square root of the power 
spectrum or "the magnitude" \N(q, u)\ of the Fourier transform of the STM measurement 
of the local conductance map. Here we will follow this convention. For the case of weak 
Born scattering, the magnitude of N(q,u) is proportional to \ImA(q,uj)/7c\. For the 2D 
electron gas, in Fig. 2, we have plotted \N(q,ui)\ normalized to y/fiiN(0)\5e\/ /jl versus q for 
various values of u/fi. Here, with q in units of 2kp one has 

w ,,„)l = , * , 1 (19) 



for (1 + w/fi) < q/2kp. As shown in Fig. 2, for the weak scattering case \N(q,Lv)\ has a 
one-sided square root singularity at q/2kp = (1 + co/n) 1 ^ 2 . This singularity is cut off and 
the response peak varies as (ikp) 1 ^ 2 when the quasi-particle mean-free path is taken into 
account. In this case 

i f r i 

\N(q,u)\ = / de = (20) 



with T = (Ikp^ 1 - Plots of \N(q,u) \ versus q for u = and several different values of T are 
shown in Fig. 3. 

Similarly, for an anisotropic system with 

k 2 k 2 

^ = 7^ + 7^-^ (21) 
2m x 2m y 

one finds that 



n (0 for q 2 x + ~f 2 q 2 < 8m x (fi + u) 



(22) 
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with 7 2 = m x /m y . In this case the cusp in the power spectrum follows a locus determined 
by 

q 2 x + l 2 q 2 y = 8m x (fi + u) (23) 

which reflects the elliptical Fermi surface. In Fig. 4, we have plotted \N(q, u) \ for to = and 
7 = 3 normalized to y/n'i{m x /2'K)\8e\/ //, which gives 

\N(q x ,q,,u)\ = ±- , , 1 , ., , 1 (24) 



with q x and q y measured in units of 2kp = 2(2m x fJ,) 1 ' 2 . Here, we see that \N(q, 0)| vanishes 
inside the ellipse Eq. (J2HJ) and has a square root divergence as q approaches the ellipse. There 
is a reduction 7 = (m^/m^.) 1 / 2 in the strength of the cusp along the q y direction relative to 
the q x direction that reflects the fact that the joint density of states which enters A (q, u) 
depends on the curvature of the Fermi surface. 

Turning now to the case of a tight-binding band, we consider first the simple nearest- 
neighbor hopping band for a square lattice with a unit lattice spacing 

efc = — 2t (cosk x + cos k y ) — \i. (25) 

The Fermi surface is shown in Fig. 5 for \i = —t, corresponding to a small filling (n) ~ 0.31 
which we have chosen to illustrate what happens when ImA(q,w) is folded back into the 
first Brillouin zone. In the following tight-binding bandstructure calculations, we will assume 
that N(x, to) is measured on a grid of points corresponding to the sites of the lattice. In this 

case, 

k 

with 

G (k) = n (27) 

00 — 6k + to sgnyuj) 

Here, N = L x L is the number of lattice sites and k and q are defined in the first Brillouin 
zone with components running from — 7r to n. We have set the lattice spacing a — 1. 
This choice of the grid simplifies the calculations but has the consequence that all results 
are folded back into the first Brillouin zone. It is this down folding that leads to the 
kaleidoscopic patterns in the numerical results shown in Ref. I8J. In the experimental FT- 

n n 

STM measurements pj, |3J, a smaller grid spacing was used leading to a g-space power 
spectrum which looks more like the extended zone picture for the q- vectors of interest. 
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Carrying out the momentum sum in Eq. (|26|). we find the results for \ImA(q,Lj)/ir\ 
shown in Fig. 6. Here in Fig. 6(a), q varies along the diagonal a-cut shown in Fig. 5 
with q x = q y . In this case, \N(q,u)\ exhibits a similar cusp to that of the free electron 
system when \q x \ > 2k X p = 2cos~ 1 (— (/i + to)/ At). This same type of behavior is shown 
as the dashed curve for the b-cut with q y = in Fig. 6(b), where we have displaced the 
numerical results by ±2tt corresponding to an extended zone scheme. Here, the cusp occurs 
for \Qx\ > 2kp = 2cos~ 1 (— (2t + fi + u)/2t) which is greater than tt for the b-cut of Fig. 5. 
In practice, the numerical data is obtained in the "reduced" (— tt, tt) zone so that the q y = 
b-cut appears as the solid curve shown in Fig. 6(b). As uo is increased, the characteristic q 
values can move across the boundary of the first zone and be mapped back via a reciprocal 
lattice vector. This can then lead to a situation in which the square root singularity is 
approached from smaller q- values with \N(q,u)\ vanishing when q exceeds a critical value 
as shown by the solid curve in Fig. 6(b). As noted, it is this folding of the FT-STM power 
spectrum into the reduced zone that leads to the kaleidoscopic 2D (q x ,q y ) patterns for 
different ui values which have been reported 

Finally, consider the case of a tight-binding band, like that found for BISCO. Here, one 
has a next-near-neighbor hopping t' so that 

€k = — 2t(cos k x + cos k y ) — 4£ cos k x cos k y — \i . (28) 

The Fermi surface for t'/t = —0.3 and fi/t = — 1.0 is shown in Fig. 7. Results of a numerical 
calculation for ImA (q,uo) /tt for this case are plotted in Figs. 8(a) and 8(b) for q y = and 
Qx — Qy, respectively. For q y = 0, the nesting vector 3 is shown in Fig. 7. If we consider 
the closed Fermi surface around (tt, 7r), one can see that the two points of this closed Fermi 
surface that are connected by the wave vector labeled 3, have Aq greater than tt. In fact, 
Aq/n ~ 1.76, so that when this is folded back into the first Brillouin zone, the cusp occurs 
at Aq/n ~ —.24. Also, as we saw previously in Fig. 6(a) the cusp rises from smaller values 
of q in the reduced zone although in an extended zone scheme it would be approached from 
larger q values just as in the free electron case. The peaks in the response for q x = q y shown 
in Fig. 8(b) arise from the nesting vectors labeled 1 and 2 in Fig. 7. 

To conclude this Section on the normal state, we consider an electron-phonon system with 
a self-energy S(c<j) which depends only upon the frequency. In this case, the momentum that 
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enters the Hankel function giving G(x, u) is 

k(u>) ^k F + {UJ ~ E(a;) . (29) 
v F 

The real part of the self-energy leads to a shift in the wave length of the modulations and 
the imaginary part of X leads to their exponential decay on a scale i = vf/{—2T, 2 (lj)). 
This behavior is reflected in a shift in the position and a rounding of the cusp in the power 
spectrum P(q, u). For values of u which are small compared to a typical phonon frequency 

k(u) ^k F + (30) 
Vf 

with A = 2iV(0)|g| 2 /f2o the dimensionless electron phonon interaction strength. Here g is the 
effective electron-phonon coupling and Qq a typical phonon energy. As u increases 2kp{oj) 
will reflect the detailed dependence of both the real and imaginary parts of £(a>). 



III. THE FT-STM POWER SPECTRUM OF A SUPERCONDUCTOR 



Next consider the case of a superconductor. Here, for an on-site energy perturbation, we 
have (4 1 

5N(x, lo) — Im 1 d 2 x' (G (x — x', u) G (x' — x, u) — F (x — x' , u) F (x' — x, u)) 5e fx') 

7T J 

(31) 



uj + e k 



with G(x,u) the usual single particle propagator 
and F(x, u) the anomalous Gor'kov propagator 



In this case jll | 



F(x,u) 



d 2 k 



d 2 k 



ik-x 



A, 



(27T 



I 2 u 2 -el-A 2 k + i5' 



(2nf 



G(k + q, u) G(k, u)-F(k + q, u)F (k, u) 



(32) 



(33) 



(34) 



For a 2D s-wave superconductor with = A and = k 2 /2m — fi, these Green's functions 
have characteristic wave vectors 



k±(u) = hi 



1 ± 



y/uj 2 - A 2 " 



(35) 
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The wave length of the ripples in SN(x,u) are set by 2k±(u). For uj < A, the ripples 
exponentially decay. If the impurity perturbation involved a change in the size of the 
gap, there would be a quasi-particle interference contribution involving k + (uj) — k_(u) = 
2kp\/uj 2 — A 2 / ji. The ripples associated with this contribution vary on a length scale set 
by Vp/A which is of order the coherence length rather than k F l . It is this type of spatial 
oscillation that is responsible for the Tomasch oscillations 12j. While Tomasch ripples at 
Qt(&) — k + (oj) — k_(uj) are not present for the case of a charge impurity in an isotropic 
s-wave superconductor, they can appear for a <i-wave superconductor as discussed below. 

As shown in Appendix B, it is straightforward to evaluate ImA (q, uj) for an s-wave BCS 
superconductor and one finds that 



1 r w x , „ r2 ^ {lj-VuP=5F\ 1 



-ImA(q,u) = 4vriV'(0) / , (36) 



for 2k^{oj) < q < 2k + {uj) and 



-ImA(q,u;) = 4vriV 2 (0) fu ~ ^ ~ ^ 



7T V Vuj 2 - A 
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+ Vuj 2 - A 2 \ 1 



1 




-4k 2 _{uj) 




1 




-4A; 2 H 




1 


Vv 2 


- 4A; 2 (uj) 



for 2A; + (u;) < g. More generally, for an electron-phonon system the frequency-dependent 
complex gap A (a;) would replace A. Normalizing by y/niN(0) \8e\/ [/,, as previously done 
for the free electron case, we have plotted \N(q,u)\ versus q in Fig. 9 for a BCS s-wave 
superconductor with A//x = 0.1. The coherence factors in Eq. ()36|) and ([37)1 give more 
weight to the 2k + {u) cusp for positive values of u associated with a bias voltage that probes 
the empty states, Fig. 9(a), while the q = 2k-(cu) cusp is enhanced when the bias is reversed 
as shown in Fig. 9(b). 

We turn next to the case of a d-wave superconductor with 

A& = A (cos k x — cos k y )/2 (38) 

and a bandstructure given by Eq. (J28|) with t'/t = —0.3 and \x = —1.0. The Fermi surface 
for these parameters is shown as a dashed line in Fig. 10(a). The contours of the solid 
regions correspond to the loci of points where u = \/e 2 + A 2 , for u = 0.5A in Fig. 10(a) 
and uj = 1.1 Ao in Fig. 10(b). Results for \ImA (q, u))/tt\ at various values of u less than Ao 
are shown in Fig. 11(a) for q x = q y and Fig. 11(b) for q y = 0. For the diagonal cut with 

11 



Qx = Qy, there is a response at the wave vector qj(uj) shown in Fig. 10(a) which connects 
the end points of the uj 2 = e\ + A\ contour. This peak is similar to the peak we have 
seen in the case of the ellipse discussed in the previous section. There are two identical 
contributions coming from the two contours on opposite sides of the Fermi surface. As the 
bias voltage eV = uj increases, one sees that ^(cu) peak moves to larger values of momentum. 
This reflects the increase in magnitude of the wavevector qj(uj) as uj is increased and can 
provide information on the fc-dependance of the Fermi surface and gap 

ELS 

. In addition, 

there is a contribution coming from c^(u}) which connects the tips of two opposite contours 
which increases more slowly with uj. Finally, there is a response associated with q' shown in 
Fig. 10(a). This is the response due to the wave vector labeled 1 in Fig. 7 and shown in the 
diagonal q x = q y response of the normal metal with this same bandstructure in Fig. 8(b). 
Approximate analytic results for the <i-wave case are given in Appendix B, Eqs. (jB14j) and 

Similarly in Fig. 11(b), for q y = one finds structure associated with the wave vectors qi 
and qs shown in Fig. 10(a). As discussed by McElroy et. al, |3j this structure arises from 
a peak in the joint density of states associated with the overlap of the ends of two opposite 
uj 2 = e 2 k + Al contours. Just as for the case of the elliptical Fermi surface previously discussed, 
the strength of the cusp depends upon the curvature of the dispersion. The increase of qi(u) 
with increasing uj again provides information on the Fermi surface and The peak at 
qs(uj) initially increases with uj and would continue to increase in an extended zone but here, 
for uj > 0.25, it is reflected back into the first Brillouin zone. Finally, there is a weak cusp 
at low momentum Qt{uj) associated with the Tomasch interference process. 

Note that the response seen in \N(q,u)\ associated with quasi-particle interference is 
actually characterized by a continuous curving cusp in the (q x ,qy) plane whose intensity is 
related to the joint density of the k and k + q quasi-particle states. When we take into 
account all four quadrants, the structure in the power spectrum P(q,u) in the first Brillonin 
zone becomes much richer. Not only are there cusps associated with scattering processes 
confined to similar contour regions in quadrants 2 and 4, which give rise to a cusp whose 
major axis is perpendicular to that due to quadrants 1 and 3, there are cusps from scattering 
processes between the four nodal regions. In addition, there is the response associated with 
the nodal wave vector q'. In order to see this more clearly, we have calculated A (q,ou) for 
the case in which t' = 0. In Fig. 12 we compare the results for \Im A (q,uj) /tt\ calculated 
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with t' — [Fig. 12(a)] and with t'/t = —0.3 [Fig. 12(b)]. In both cases, we clearly see 
the contributions of q3 and q 7 . In addition, a q' contribution from the nodal region is also 
visible. The q 3 contribution is weaker for the case in which t' = 0. In this case we also find 
that | q 3 1 < |q'| as one can see from the insets in the figures. 

When uj is greater than A , structure in \N(q, uj)\ is associated with the Q a vectors shown 
in Fig. 10(b), for uj = 1.1A . In this case, for q x = q y , \N(q, uj)\ exhibits structure at q x = Qi, 
Q2, and Q3, see Fig. 13(a). The weak structure at large momentum transfer Qi is similar 
to the structure for the normal state labeled "1" in Fig. 8, since the gap vanishes along the 
diagonal. The structures at Q2 and Q3 correspond to interference processes between the 
particle like and the hole like BCS quasi particles. Just as for the s-wave case, for uj > 0, 
the coherence factor gives more weight to the Q 2 process. Results for \N(q, cu)\ versus q x 
for q y = are shown in Fig. 13(b). In this case, structure appears associated with q x = Q4 
and Q 5 with the coherence factors leading to a large response at Q4 for positive uj. As 
noted earlier, in the cuprate superconductors lifetime effects suppress this uj > A structure, 
making it difficult to see in the experimental studies. 

The structure in the quasi-particle interference response can also be seen in the intensity 
map plots of \ImA (q,uj)/ir\ over the (q x ,q y ) plane. One such map for uj = 0.5A is shown 
in Fig. 14. Here one sees that the response is characterized by continuous curving intensity 
cusps in the (q x , q y ) plane. Like the case of the elliptical Fermi surface, the intensity can have 
significant variations along the cusps due to the curvature of the quasi-particles dispersion 
relation. Going out from the origin along q x = q y as in Fig. 12(b), one first sees a bright 
(high intensity region) cusp associated with q 7 . At large momentum values along this same 
45° line one sees a narrow bright line associated with q' scattering processes which connect 
the outter edges of the uj = ^ e\ + A| envelopes. Finally, the brighter, curved region of 
intensity near the (tt, it) corner arises from q 3 interference processes associated with the 
inner boundaries of these contours. The bright regions near (q x /7i = 0.85, q y /ir = 0.3) 
and (q x /n = 0.3, q y j V = 0.85) arise from interference effects associated with g 2 and q% of 
Fig. 10(a). 

Finally, we consider the response of a <i-wave superconductor when there is an ordered 
stripe array of weak scattering centers. In Fig. 15 we show the ripples produced in N(x, uj) 
when the scattering centers form a stripe-like structure oriented along the y-axis with a 
spacing of four lattice sites. For the parameters we have chosen, q^u) ~ n/2 for uj = 0.5A 
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with A the maximum value of the <i-wave gap. As seen in Fig. 15, the ordered array of 
stripe scattering centers produces an oriented set of ripples. These give rise to the structure 
in \N(q, u)\ shown in Fig. 16 for q ~ n/2. Here we have assumed that there is quasi-particle 
scattering from a random array of impurities as well as the stripes. These contributions add 
incoherently so that 



<l^)| 2 > 



2\ 



1 + A- 



(39) 



~(q x -Q x ) 2 + T\ 

Here Q x = n/2 and we have taken T = 0.01 and A = 0.5 in Fig. 16 to illustrate the possible 
interplay of the scattering from the random impurities and the stripes. The quasi-particle 
interference peak associated with qi moves to lower values of q x as u increases while the 
response associated with the ordered array of scattering centers remains fixed at q x ~ 0.57T 
and only its amplitude changes with uj. For larger values of A, the stripes would be the 
dominant feature while for smaller values of A scattering from the random impurities would 
dominate the power spectrum. 



IV. CONCLUSIONS 



Here we have discussed some detailed examples which illustrate what can be learned 
from FT-STM studies of layered materials. For the case of a weak potential perturbation 
we have seen that the wave vector power spectrum of the local tunneling density of states, 
P(q,u), contains information on the quasi-particle spectrum and the structure factor of the 
scatterers. This is also the case for dilute impurities even if they act as strong scattering 
centers and Se is replaced by a t-matrix. For a normal metal, one can determine information 
about the nesting properties of the Fermi surface from the loci of the q-space cusps for 
uj — > 0. In addition, the dressed Fermi velocity can be obtained from the u dependence 
of the position of the cusp. In fact, from the uj dependence of the cusp position and its 
rounding, one can, in principle, obtain information on the real and imaginary parts of the 
electron self-energy. 

In the superconducting state, one can also obtain information on both the wave vector 
and frequency dependence of the gap A(k, uj). Again, just as in the case of a normal metal, 
the signature of these quasi-particle interference effects are continuous "arcs" of cusps in 
q-space. However, the intensity variation along the arcs may be large due to changes in the 
effective density of states associated with a given momentum transfer q, making the intensity 
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appear more like spots. The rapid change in intensity is due to the small parameter Ao/t. 

In addition to these effects, one can also obtain information on the static structure factor 
of the scattering potential (\5e q \ 2 ). As discussed in Appendix A, for the case of a local 
potential, averaging over a Aq x by Aq y block of q values about a given q leads to (\5e q \ 2 )/N ~ 
riiSe 2 . However, if the scattering potential has long-range order, such as is the case for the 
B1O2 over layer in the BISCO studies, one should see a g-dependent response in the wave- 
vector power spectrum of the local tunneling density of states for q values equal to the 
reciprocal lattice vectors of the B1O2 layer. This response will be particularly strong for 



uj ~ Aq. Similarly, it shou 



D 



d be possible to see evidence of pinned stripes in the structure 



factor if they are present 

As discussed in the introduction, we have focused on the case of weak scattering. As 
noted, in this limit, where the Born approximation is adequate, one has a simple separation 
of the response into the quasi-particle interference effects which are contained in A(q,u) and 
the structure factor of the scatterers which is contained in (|5e g | 2 ). Of course, many-body 
interactions will give rise to additional effects. There will be screening, changing 5e q to 
5e q /e(q, 0) where e(q, 0) is the zero frequency dielectric constant j^]. Furthermore, there will 
be vertex corrections so that for example, for a strongly-interacting normal system one will 
have 

A(g, u)=J -0^ T (k, g, u) G(k + q, u) G(k, u) , (40) 

with F(k,q, to) the elastic vertex for momentum transfer q and zero energy transfer for an 
incoming particle with momentum k and energy uj. Similar vertex corrections will occur in 
the superconducting state. In addition, the quasi-particle dispersion relation can be altered 
by an interaction which breaks the translational symmetry leading to a different uj 
dependence of \N(q, uj) \. 

Finally, there is the form factor of the tunneling probe. Here we have neglected the 
momentum dependence of the tunneling matrix element and simply assumed that the con- 
ductance map was proportional to the local tunneling density of states N(x,u). However, 
if this is not the tunneling matrix element form factor 

T(fc) = ^V fc -'T(/) (41) 
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will enter so that for a superconductor 

For example, a tunneling form factor 

T(fc) = (cos /c x — cos k y ) (43) 

n 

has been suggested for the case of an STM tip on a cuprate superconductor p. However, 
because the average tunnelling density of states in these STM experiments appears to vary 
linearly with V at low voltages, we have modeled the tunnelling as a direct process and 
neglected its momentum dependence. 



APPENDIX A: RANDOM IMPURITY STRUCTURE FACTOR 

Suppose we have an L x L section of a lattice with a concentration of impurities per 
unit area. Assume that if there is an impurity at site I it has a potential 5e(l) = 1, while if 
there is no impurity at site I, 5e(l) = 0. For a given random configuration of iVj impurities 
on L x L = N sites, corresponding to an area impurity concentration n« = Ni/N, we have 

fe( g ) = $>*'fc(0- (Al) 
i 

For this configuration of impurities we can define a power spectrum with 

P{3 ) = ^f- (A2) 
There is of course a peak for q = 0, where 

P(0) = -f- (A3) 



However for other values of q, P{q) fluctuates about n,. 
many realizations of independent impurity configurations 



i one were to average P(q) over 

3, 



(P(q)) = m for q ^ 0. (A4) 

However, if we have an L x L lattice with a fixed configuration of impurities, the power 
spectrum P(q), given by Eq. (jA2|) . will exhibit fluctuations. If we average over blocks 
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of q of width (Aq x ,Aq y ) about each q we can reduce these fluctuations. Naturally at the 
same time, the momentum resolution will be reduced. Define a block smoothed P(q) average 
as follows: 



here the prime on the sum indicates that -P(O) is replaced with rij so that the q = peak is 
not broadened from the smoothing operation. The sum in Eq. (jA5|) is over a set of q-points 
(—Aq x /2, Aq x /2) and (—Aq y /2, Aq y /2) about each q. N(Aq x , Aq y ) is the number of sites in 
the Aq x by Aq y block. We expect that the RMS fluctuations of P(q) will decrease inversely 
as the square root of the number of g-values in the Aq x by Aq y block. 

For example, consider a lattice with L = 500 and one configuration of random impurities 
with rii = 0.01. The power spectrum P(q), given by Eq. (|A2j) . versus q x for q y = is shown 
in Fig. 17(a). Figures 17(b) and 17(c) show similar plots for the smoothed power spectrum 
P(q), given by Eq. (|A5|) . for square blocks with Aq x = Aq y = 2ir{2d + l)/L for d = 5 and 
d = 10 respectively. Figure 17(d) shows P(q) for the case in which only Aq y = 2ir(2d+ 1)/L 
is averaged with d = 10. The solid line in each of these figures is rij = 0.01. 

It is clear from the results shown in Fig. 17, that one can reduce the fluctuations in the 
power spectrum by averaging it over a region of q-space surrounding a given q-point. This 
is of course just complementary to taking the original L x L spatial lattice, breaking it up 
into Lb x Lb blocks with Lb = 2ir/Aq, constructing the appropriate Fourier transform for 
each block and then averaging over the blocks. The more blocks one has, the small the RMS 
fluctuations of the power spectrum of the Fourier transform. Of course, dividing the L x L 
system up into more blocks leads to a corresponding decrease in the momentum resolution. 
The momentum smoothing operation, Eq. (|A5j) . is just another way of blocking. 

As a further test of these ideas, we calculate the RMS deviation of P(q) from versus 
the momentum block size 2d + 1 



In Fig. 18 we show that D varies as (2d+ l) -1 for large d for several different concentrations 
of impurities. The straight lines on the plot represent the asymptotic behavior in which 
D ~ arii/ {2d + 1) with a of order one. 




(A5) 




(A6) 
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Finally, we model the behavior of the magnitude of the wave-vector Fourier transform of 
the local tunneling density of states by 

if q > q* = it/A 



\N(q)\ = y/Pfijx 1 (A7) 

I otherwise. 

The results for L = 500 and rij = 0.01 are shown in Figs. 19(b)-(c) for the 
BlockAvg [\N(q x: q y = 0)|], averaged over a Aq x x Aq y square with d — 5 and d — 10 
respectively. Figure 19(a) shows the unaveraged result. Similarly Fig. 19(d) shows the 
LineAvg [\N(q x , q y = 0)|] versus q x for the case in which we only average over a q y segment 
of width Aq y = 2n(2d + 1)/L with d = 10. Thus the g-averaging brings out the underlying 
structure of the quasi-particle interference factor. 

Finally, we note that even without such block averaging, one is still able to see an image 
of the quasi-particle interference response in an intensity plot of 



\N(q,u>)\ = 



—Ira A (q, u) 

7T 



Mq)l (A8) 



'N 

over the (q x ,qy) plane. In Fig. 20 we show a plot of \N(q, u)\ for the same parameters 
that were used in Fig. 14. Here |<5e(g)| is obtained from Eq. (Al) with one realization of 
an impurity concentration rii = 0.01. Comparing this figure with Fig. 14, one can see that 
while the structure is noisy, the key features remain clearly visible. 

APPENDIX B: QUASI-PARTICLE INTERFERENCE RESPONSE 

The quasi-particle interference response function for the 2D free electron gas 
1 1 r 

-ImK{q,u) = -Im / d 2 x e* qx G (x, u) G a (-x, u) (Bl) 

7T IX J 

can be directly evaluated using the expression for G (x,cu) given by Eq. (15). Since the 
Hankel function = J + i N , we have 

^-ImA(q,uj) = -2nN\0) J d 2 xe iq ' x J Q (k(cu)r) N (k{u)r) 

8ttA^ 2 (0) i . 1 = q>2k(u) 



with r = \x\ and 



q < 2k(uj) 



h-U) =kr i l + J L ) ( B3 ) 
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We also note that Carrying out the angular integration one has 

q > 2k(u) 



8nN(0) 2 



ReA(q, to) 



q 2 y/q 2 -4k 2 



hi 



y/q 2 -4k 2 -l 



-167riV(0) 2 



tan" 1 9 q < 2k(u) 



(B4) 



Plots of the real and imaginary parts of A(qu) are shown in Fig. |^ As noted, in Born 
approximation one only sees Im\(q,oj). However, when the scattering is stronger, giving 
rise to a phase shift, both Re\(q,u>) and Im\(q,u) will be present in the quasi-particle 
interference response. The imaginary part of A can also be evaluated by noting that for 

LU > 0, 

d 2 k 1 1 



1 1 f 

— ImA(q,uj) = — Im I 
IX ix J 



(27r) 2 e k+q - e k 

Carrying out the angular integration one has 



(27r) 2 UJ — t k+q + l5 CO — Cfc + l5 

,H ' 5{u-e k ) (B5) 



2- 



2rn 



q 2 > Ak 



y/q 4 ~4k 2 q 2 /ggN 



so that 



2vr e k+q - e k } q q 2 < Ak 2 

q 2 /8m— fj, 

1 f 2m 

-ImA(q,uj) = 2N(0) / de k 6(u - e k ) - = (B7) 
TT J Jq A - Ak 2 q 2 

which gives Eq. ()B2j) since N(0) = m/2ix. 

For the electron-phonon case in which the self-energy £(a>) = +2X2(0;) only depends 

on u, the angular integral can be evaluated in the same way leading to 

q 2 /8m~[i 

-ImA(q,uj) = 2N(0) f de k 2m (B8) 

Here, as usual for the electron-phonon problem we have neglected vertex corrections. 
For the case of a superconductor with scattering from a site charge potential 

A ^) = / (27F c^-^ + i*) p^TTs) (B9) 

with = -^/e| + A|. For a constant s-wave gap A^ = A and 

1 r a / x 2 T [ d 2 k uj + e k 1 
ImA g,u, =— / -— ^ 2 " BIO 
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7T 7T 



g 2 /8m-/i 



Making use of Eq. (|B6|) , we have 

1 f 2m 

-ImA(q,u) = 2N(0) / de k 8 (u 2 - E 2 k ) (u + e k ). (Bll) 

Carrying out the e k integration leads to the results given by Eqs. - (|S7jl in the text. 
For a ci-wave superconductor, we have 



- j m A M = 2 y ^ f („ - fr) — ^- ^ (B12) 

This integral can be approximately evaluated when co> < A for the case in which e k = 
k 2 /2m — k 2 F /2m and A^ = Aocos(2#). For example, for q along the 45° direction with 
q = Qx = % and q > kpUj/Ao, one finds that 

— 7mA (g, cj) ~ iV 2 (0) ( - ) - L (B13) 



VAT(0)A o ; ^q-^ 
This corresponds to the contribution which comes from a momentum transfer that connects 



the ends of one u = yj e 2 k + A 2 contour (i.e. a qy(tu) wave vector) similar to the case of the 
ellipse discussed in Section III. The enhancement factor (iV(O)Ao) -1 arises from the large 
curvature and resulting large density of states at the contours ends. 

Keeping q along the 45° diagonal where q = q x = q y , there are additional square root 
peaks in Im A(q, w)/tt which arise when q connects two different contours. There is a q3 (un- 
like peak near which 

1 T a/ x 4vriV 2 (0) / u \ 1 1 

-ImA(q,u) ~ — — — - (B14) 

1 ' q V2A ; Vki^q y/^u) - q 1 ' 

when ^ 

q<q z {u) = 2k F \l-\ I. (BIO) 



2 V 2A 



There is also a peak associated with 



2 



where 



q'(uj) = 2k F + (B16) 

1 , . . 47riV 2 (0) f u \ 1 1 

-ImA(q,u)~ — — — ; B17 

q V2A ; y/ q - q'(uj) v/g-g 3 M 



when q > q'{oj). In the limit Aq — »■ 0, this last expression becomes 



1 4ti-N 2 (0) 

-ImA(q,u)~ W (B18) 

n qyk F \/q - 2k[uj) 



which is the free electron result Eq. flB2|) for u/fi <C 1. 
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FIG. 1: (a) The local tunneling density of states N(x,u) due to scattering from a weak potential 
at the origin y], for (a) u = 0.5Ao and (c) u = l.lAo- The length scale is set by kp 1 . These results 
were obtained for a cylindrical Fermi surface with a d-wave gap A(#) = Aqcos(2#). The contours 



of the solid regions in (b) and (d) correspond to the points where uj = y e 2 . + A 2 (6) for uj = 0.5Ao 
and l.lAo, respectively. The dashed line in (b) is the Fermi surface for the non-interacting system. 
The q a (w) wave vectors in (b), introduced in Ref. connect the tips of various contours and 
q'(u;) connects what would be a diagonal nesting vector for the normal state Fermi surface. The 
Q a (cj) wave vectors shown in (d) are the relevant nesting vectors when uj > Aq. 
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FIG. 2: Plot of the normalized \N(q,u>)\ versus q for a free electron gas for various values of oj//j,. 
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FIG. 3: Plot of the normalized mean-free path corrected \N(q,u)\ versus q for an electron gas with 
lo = 0. Here V = (ikp)- 1 : V = (solid), T = 0.01 (long dashed), T = 0.05 (dashed), T = 0.1 
(dotted). 
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FIG. 4: Plot of the normalized \N(q,oj)\ versus q, Eq. (|24|) . for an elliptical Fermi surface with 
7 2 = m x /m y = 9 and uj = 0. In (a) q x = and in (b) q y = 
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FIG. 6: \Im A(q,uj)/ir\ for the tight-binding band, Eq. l25| with t = 1 versus q for u = and (a) the 
diagonal cut-a and (b) the horizontal cut-b indicated in Fig. 5. The dashed curves show the results 
in the second Brillouin zone. For the cut-a, 2k F y ~ .8ir so the solid curve in the first Brillouin looks 
similar to the free electron response. However, for cut-b, 2kp ~ 1.37T so that the dashed curve in 
(b) looks similar to the free electron response and the solid curve in the first Brilloin zone results 
from folding the dashed curve back into the first zone. 
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FIG. 8: \Im A(q, oj)/it\ versus q for oj = for the tight-binding Fermi surface shown in Fig. 7. (a) 
q y = and (b) g x = q y . 
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FIG. 9: Plot of the normalized \N(q,u)\ versus q for an s-wave superconductor with a cylindrical 
Fermi surface and A/fi = 0.5. (a) Results for positive values of lj = eV and (b) for negative values 
of UJ. 
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FIG. 10: Characteristic wave vectors for a d-wave superconductor with a tight-binding band, 
Eq. (|28|). with i'/i = —0.3, fj,/t = —1.0, and Ao = O.li. The contours of the solid regions 
correspond to the points where to = ye^ + with (a) u = 0.5Ao, and (b) u> = l.lAo- The 
dashed line in (a) shows the Fermi surface of the non-interacting system. In (b) Q3 and Q§ (Q4 
and Q2) connect the innermost (outermost) surfaces while Q\ connects the outermost surfaces. 
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FIG. 11: \Im A(g,o;)/7r| versus q for a d-wave superconductor with a tight-binding band given by 
Eq. ()28|) with i'/f = — 0.3, /i/i = —1.0, Ao = O.li. Here (a) shows results for a diagonal cut for 
Qx = <7y and (b) shows results for a horizontal q y = cut. The g Q and g' vectors are shown in 
Fig da). 
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FIG. 12: \Im A(q, uj)/it\ versus q for a diagonal cut in which q x = q y . (a) shows results for a 
tight-binding band with t' /t = and (b) shows results for t' /t = —0.3. to = 0.5Ao and the other 
parameters are the same of Fig. 1111 
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FIG. 13: \ImA(q,uj)/ir\ versus q for a ci-wave superconductor with a tight-binding band given by 
Eq. (gSJ with t'/t = -0.3, fi/t = -1.0, A = O.lt, and a; = 1.1A . Here (a) shows results for a 
diagonal 45° cut for q x = q y and (b) shows results for a horizontal q y = cut. The Q a vectors are 
shown in Fig llOf b). 
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FIG. 14: Intensity map of the quasi-particle interference response \Im A(g, cj)/tv\ for a d-wave 
superconductor plotted over the first (q x ,Qy) quadrant. Here, the band structure is given by 
Eq. (j2EJ) with t'/t = 0.3 and ji = -1.0, A k is given by Eq. (38) and the bias is u = 0.5A . 
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FIG. 16: Structure in |iV(g,w)| versus q x for random impurities and an ordered array of stripes 
separated by 4 lattice spacings. The quasi-particle interference peaks associated with qi(u) and 
95(0;), previously shown in Fig. 11(b), are seen along with a peak at q ~ 0.5-7T which arises from 
the striped array of scattering centers. 
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FIG. 17: Power spectrum, as defined in Eq. (|A2|) for a random configuration of impurities with 
concentration rii = 0.01. Here (a) shows the unaveraged results, (b) and (c) show the average over 
blocks of size Aq x = Aq y = 2ir(2d + 1)/L with d = 5 and d = 10 respectively, and (d) shows the 
average over lines of length Aq y = 2ir(2d + 1) with d = 10. 
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FIG. 18: RMS deviation of P(q) from m as defined in Eq. ()A6|) . The straight lines represent the 
asymptotic behavior in which D ~ otrii/(2d + 1). 
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FIG. 19: The magnitude of the Fourier transform of the local tunneling density of states, as defined 
in Eq. (|A7|) for a random configuration of impurities with concentration ni = 0.01. Here (a) shows 
the unaveraged results, (b) and (c) show the average over blocks of size Aq x = Aq y = 2ir(2d+l)/L 
with d = 5 and d = 10 respectively, and (d) shows the average over lines of length Aq y = 2n(2d+ 1) 
with d = 10. 
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FIG. 21: The real (solid) and imaginary (dashed) part of the quasi-particle interference response 
versus q for a non-interacting 2D electron gas. 
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